% 2-D BIBLIS Test Problem, keff =  1.025110, I get 1.025124 w/ extended mod

it = 1;

dat = [ 1.4360000  0.0272582  0.0058708  1.0000000  0.0000000  0.0000000
        0.3635000  0.0750580  0.0960670  0.0000000  0.0177540  0.0000000
        1.4366000  0.0272995  0.0061908  1.0000000  0.0000000  0.0000000
        0.3636000  0.0784360  0.1035800  0.0000000  0.0176210  0.0000000
        1.3200000  0.0257622  0.0000000  0.0000000  0.0000000  0.0000000
        0.2772000  0.0715960  0.0000000  0.0000000  0.0231060  0.0000000
        1.4389000  0.0274640  0.0074527  1.0000000  0.0000000  0.0000000
        0.3638000  0.0914080  0.1323600  0.0000000  0.0171010  0.0000000
        1.4381000  0.0272930  0.0061908  1.0000000  0.0000000  0.0000000
        0.3665000  0.0848280  0.1035800  0.0000000  0.0172900  0.0000000
        1.4385000  0.0273240  0.0064285  1.0000000  0.0000000  0.0000000
        0.3665000  0.0873140  0.1091100  0.0000000  0.0171920  0.0000000
        1.4389000  0.0272900  0.0061908  1.0000000  0.0000000  0.0000000
        0.3679000  0.0880240  0.1035800  0.0000000  0.0171250  0.0000000
        1.4393000  0.0273210  0.0064285  1.0000000  0.0000000  0.0000000
        0.3680000  0.0905100  0.1091100  0.0000000  0.0170270  0.0000000];

numg = 2;  numm = 8;
w=23.1226; % assembly width in cm
xcm  = [  0  0.5*w  1.5*w  2.5*w  3.5*w  4.5*w  5.5*w  6.5*w  7.5*w  8.5*w];
xfm  = [  10    20    20    20    20    20    20    20    20    ]/1;
ycm  = [  0  0.5*w  1.5*w  2.5*w  3.5*w  4.5*w  5.5*w  6.5*w  7.5*w  8.5*w];
yfm  = [  10    20    20    20    20    20    20    20    20    ]/1;
    
mt  = [  1  8  2  6  1  7  1  4  3  
         8  1  8  2  8  1  1  4  3
         2  8  1  8  2  7  1  4  3
         6  2  8  2  8  1  8  4  3 
         1  8  2  8  2  5  4  3  3
         7  1  7  1  5  4  4  3  3
         1  1  1  8  4  4  3  3  3
         4  4  4  4  3  3  3  3  3
         3  3  3  3  3  3  3  3  3 ];    
BCL = 1;
BCR = 0;
BCB = 1;
BCT = 0;      

src(:,:,1) = 0;        
src(:,:,2) = 0;        
%src(:,:,3) = 0; 

IBSL = [0 0 0 0 0
      0 0 0 0 0]';
IBSR = [0 0 0 0 0
      0 0 0 0 0]';
IBSB = [0 0 0 0 0
      0 0 0 0 0]';
IBST = [0 0 0 0 0
      0 0 0 0 0]';